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Abstract 

We report on a preliminary measurement of time-dependent CP-violating asymmetries in B — > 
J/ip Kg and B° — > t/j(2S)Kg decays recorded by the BABAR detector at the PEP-II asymmetric 
B Factory at SLAC. The data sample consisted of 9.0 fb _1 collected at the T(45) resonance and 
0.8 fb -1 off-resonance. One of the pair of neutral B mesons produced at the T(4S) was fully 
reconstructed, while the flavor of the other neutral B meson was tagged at the time of its decay. 
The value of the asymmetry amplitude, sin2/3, was determined from a maximum likelihood fit 
to the time distribution of 120 tagged candidates to be sin2/3 = 0.12 ± 0.37(stat.) ± 0.09(syst.) 
(preliminary) . 
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The CP-violating phase of the three-generation Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix can provide an elegant explanation of the CP-violating effects seen in decays of 
neutral K mesons [Q]. The unitarity relations between the elements of the CKM matrix can be 
expressed as six triangles of equal area in the complex plane. A nonzero areaQ directly implies 
the existence of a CP-violating CKM phase. The most experimentally accessible of the unitarity 
relations involves elements V u \, and Vtd, and is known as the Unitarity Triangle, where angles, and 
hence size of CP- violating asymmetries are expected to be large^]. Observing the CP-violating 
asymmetry can thus provide a crucial test of the Standard Model. 

The CP-violating asymmetry in b — ► ccs decays of the B° meson such as B° /B° — ► J/ipK® (or 
B° I B° — ► ip(2S)K®) is caused by the interference between mixed and unmixed decay amplitudes. 
With little theoretical uncertainty, the phase difference between these amplitudes is equal to twice 
the angle (3 of the Unitarity Triangle. In e + e~ storage rings operating at the T(4S) resonance a 
B°B° pair produced in T(45) decay evolves in a coherent P-wave until one of the B mesons decays. 
For this measurement, one of the B mesons (Bqp) is fully reconstructed in a CP eigenstate {J/i)K° s 
or ip(2S)K®). If the other B meson (B tag ) is determined to decay to a state of known flavor at a 
certain time tt ag , Bqp is at that time known to be of the opposite flavor. By measuring the proper 
time interval At = tcp — t tag from the B tag decay time to the decay of the Bqp, it is possible to 
determine the time evolution of the initially pure B° or B° state. The experimental time-dependent 
decay rate into the Bqp final state is given by 

T± = ire _r|At| [1 ± Vsin2p x sinAm d At] ®K(At; a) , (1) 

where the + or — sign indicates whether the Bt ag is tagged as a B° or a B°, respectively. The 
"dilution factor" T> is given by T> = l — 2w, where w is the mistag fraction, i.e., the probability that 
the flavor of the tagging B is identified incorrectly. The term 1Z accounts for the finite detector 
resolution, where a represents the set of parameters that describe the resolution function. 

1 Sample Selection 

For this analysis we used a sample of 9.8 fb -1 of data recorded by the BABAR detector || between 
January 2000 and the beginning of July 2000, of which 0.8 fb -1 was recorded 40MeV below the 
T(45) resonance (off-resonance data). At PEP-II, B mesons are produced in the asymmetric 
collisions of 9 GeV electrons and 3.1 GeV positrons, and have an average boost along z direction 
of (/3 7 ) = 0.56. 

J/ip candidates were identified through their decays into e + e~ and // + //~, and ip(2S) candi- 
dates were reconstructed in e + e~, fi + fJ>~, and J/ipir + Tr~ modesQ. In the electron modes of J/ip 
and ip{2S), at least one electron was required to be positively identified. An algorithm for the re- 
covery of Bremsstrahlung photons j7j was used if both electrons were positively identified. Similarly, 
for the candidates, at least one muon was required to pass particle identification criteria. 

The charmonium states then were selected in a window around their expected mass[[IJ, and were 
combined with K® candidates reconstructed through their decays into tt + tt~ and states 0] to 
form Bcp candidates. Bqp candidates were identified with a pair of nearly uncorrelated kinematic 
variables: the difference AE between the energy of the Bqp candidate and the beam energy in 
the center-of-mass frame, and the beam-energy substituted mass msE 0- The signal region was 
defined by 5.270 GeV/ c 2 < msE < 5.290 GeV/ c 2 and an approximately three-standard-deviation cut 
on AE (typically \AE\ < 35MeV). The CP sample used in this analysis was composed of 168 
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candidates: 121 in the J/i>K° s (K% -> vr+vr-) channel, 19 in the J/tpK® (K% —* 7r°7r°) channel and 
28 in the i/>(2S)K% (K% -► tt + tt~) channel. 

2 Time Resolution Function 

All tracks in the event, excluding those from Bqp, were fit to a common vertex using an iterative 
procedure, which at each step removes a track with the worst x 2 greater than 6. Information on 
the direction of the Bqp, derived from the kinematics of the reconstructed Bcp and the average 
position of the interaction point, was also used in the fit. In order to reduce bias, all pairs of tracks 
that can be reconstructed as long-lived neutral particles (V°s), were replaced by the corresponding 
V° candidates in the fit. Events were rejected if either Bcp or Bt ag does not converge, or if 
reconstructed Az or its uncertainty were large(|Az| > 3 mm or a^ z > 400 fj,m). 

The time resolution function is described accurately by the sum of three Gaussian distributions, 
which has six independent parameters: 

3 , 

K{At;a) = ]T — ^= exp (-(At - ^Ax* 2 ) • (2) 

A fit to the time resolution function in Monte Carlo simulated events indicated that most of 
the events (/i « 1 — ji = 70%) were in the core Gaussian, which has a width a\ 0.6 ps. The 
wide Gaussian had a width a<i ~ 1.8 ps. Tracks from forward-going charm decays included in the 
reconstruction of the Bt ag vertex introduced a small bias, 5\ ~ — 0.2ps, for the core Gaussian. The 
third Gaussian, used to parameterize events at very large values of Az due to poorly reconstructed 
vertices, had a fixed width of 8 ps and a fraction f w ~ l%f|. In likelihood fits, the widths of the 
first and second Gaussian were parameterized by event-by-event error a At from the vertex fits and 
two scale factors 5i and 52 (oi = S\ X a At and o<i = Si x a At)- Parameters S2 = 2.1 and fi = 0.25 
were determined from Monte Carlo, bias of the second Gaussian 62 was fixed at ps, and three 
remaining parameters a = {Si,6i, f w } were determined from the observed vertex distributions in 
the hadronic B° decays. The time resolution is dominated by the precision of the Bt ag vertex 
position, and we found no significant differences in the Monte Carlo simulation of the resolution 
function parameters for the various fully reconstructed decay modes ||. 

3 B flavor tagging 

Each event with a CP candidate was assigned a B° or B° tag if the rest of the event (i.e., with 
the daughter tracks of the Bqp removed) satisfies the criteria for one of several tagging categories. 
Three tagging categories relied on the presence of a fast lepton or one or more charged kaons in 
the event. Two categories, called neural network categories, were based upon the output value of a 
neural network algorithm applied to events that have not already been assigned to lepton or kaon 
tagging categories. 

The mistag fractions W{ were measured directly for each tagging category by studying the prob- 
ability of flavor mixing in events in which one B° candidate, called the B rec , was fully reconstructed 
in a flavor eigenstate modeQ. Both integrated and time-dependent analyses of the fraction of mixed 
events are sensitive to the to the mistag fraction. In the time-dependent analysis, the probability 
density functions of unmixed (+) and mixed (— ) events were defined similar to Eq. (|l|) 

H± = ire _r|At| [1 ± V x cosAmrfAt] ®K( At; a) . (3) 
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These functions were used in an maximum maximum likelihood fit to the large samples of fully 
reconstructed B° decays into hadronic and semileptonic modes to determine mistag fractions W{ = 
|(1 — T>i) for each tagging category i. We find the tagging efficiency e = (76.7 ± 0.5)% and the 
effective tagging efficiency Q = e (1 — 2w) 2 = (27.9 ± 0.5)% (statistical errors only)Q. Out of 168 
CP candidates, 120 were tagged: 70 as B° and 50 as B°. 

4 Extracting sin2/3 

The value of sin2/3 was determined from an unbinned maximum likelihood fit to the time-dependent 
distribution of tagged Bqp decays ||. The fitting procedure was extensively tested in "toy" Monte 
Carlo simulations, where effects backgrounds, different resolution function parameterizations, and 
behavior of the fit on small statistical samples were studied. We also verified the reconstruction 
and fitting procedure on a large sample of events produced with the full BABAR GEANT3 detector 
simulation. 

The analysis was also validated on the various data control samples where the CP asymmetry 
was expected to be zero. An "apparent CP asymmetry" was studied on a sample of B + — > J/ipK + 
events and events with self-tagged B° — > J/ipK*° (K*° — > K + ir~) decays. We also use the event 
samples with fully-reconstructed candidates in charged or neutral hadronic modes. In all cases, the 
apparent CP asymmetry was found to be consistent with zeroj|]. 

We have adopted a blind analysis for the extraction of sin2/3 in order to eliminate possible 
experimenter's bias. We used a technique that hides not only the result of the unbinned maximum 
likelihood fit, but also the visual CP asymmetry in the At distribution. The error on the asymmetry 
was not hidden. With these techniques, numerous systematic studies could be performed while 
keeping the numerical value of sin2/3 hidden. The analysis procedure for extracting sin2/3 was 
frozen, and the data sample fixed, prior to unblinding. 

Using the maximum-likelihood fit to the full data sample of B°/B° — > J/ipK® and B°/B° — > 
ip{2S)K° s events, we determined pi 

sin2/3 = 0.12 ± 0.37(stat.) ± 0.09(syst.) (preliminary). (4) 

Systematic errors arise from uncertainties in input parameters to the maximum likelihood fit, 
incomplete knowledge of the time resolution function, uncertainties in the mistag fractions, and 
possible limitations in the analysis procedure. 

5 Conclusions and prospects 

We have presented the first measurement of the CP-violating asymmetry parameter sin2/? in the B 
meson system by BABARCollaboration. Our measurement is consistent with the presently available 
data||, and also agrees within two standard deviations with the value determined by other con- 
straints of the Unitarity Triangle^]. While the current experimental uncertainty on sin2/3 is large, 
the next few years will bring significant improvements in precision. Measurements of CP-violating 
asymmetries in other final states are also underway, and will help constrain the Standard Model 
description of CP violation. 
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